PM/97-47 
September 1997 



; Dynamical Chiral Symmetry Breaking from 

Variationally Improved Perturbative Expansion [] 

> 

q ' J.-L. Kneur 

(N 



Physique Mathematique et Theorique, UMR-CNRS, 
Universite Montpellier II, F-34095 Montpellier Cedex 5, France. 

> 
O 

o 

CO 
^— H 

We review how a specific resummation of the so-called "delta-expansion" , applied to the 
QCD Lagrangian, transforms the ordinary perturbative expansion in a s into an expansion 

r£3 ' in an arbitrary mass parameter, around the basic scale A. When applied to the pole mass, 

the resulting expression may be interpreted as a dynamical mass ansatz, to be optimized 
with respect to the new expansion (mass) parameter. The construction is generalized to 
obtain estimates of the order parameters of the SU{nf)i x SU(nf)n (n/ = 2, 3) symmetry 

^ ' breakdown. 
> 

X 
S3 



x Talk given at the QCD97 Euroconference, July 1997, Montpellier France, to appear in the proceedings. 
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We review how a specific resummation of the so-called "delta-expansion", applied to the QCD Lagrangian, 
transforms the ordinary perturbative expansion in a 3 into an expansion in an arbitrary mass parameter, around 
the basic scale A. When applied to the pole mass, the resulting expression may be interpreted as a dynamical mass 
ansatz, to be optimized with respect to the new expansion (mass) parameter. The construction is generalized to 
obtain estimates of the order parameters of the SU(rif)L x SU(nf)n (rif = 2, 3) symmetry breakdown. 



1. Introduction 

There is a common lore that a "first principle" 
determination of the order parameters character- 
izing the (chiral) dynamical symmetry breaking 
(DSB), such as the (qq) condensate typically, is 
definitively out of the reach of the basic QCD 
perturbation theory. This is largely justified, tra- 
ditionally, by the fact that DSB is an essentially 
non-perturbative mechanism. However, it may 
depend on what exactly one means by perturba- 
tion theory. For instance, since the pioneering 
work of Nambu and Jona-Lasinio (NJL) it 
has been understood how it is possible to resum 
a relevant class of graphs to obtain the qualita- 
tive (and some quantitative as well) properties 
of DSB explicitly, at least in specific approxima- 
tions and/or models. Also, independently of the 
NJL idea, the modern Chiral Perturbation The- 
ory (ChPT) (|] gives a consistent effective de- 
scription of data at low energies where the QCD 
perturbative series is not applicable. Indeed, defi- 
nite progress have been made to relate ChPT with 
generalized NJL models ||, although a precise 
connection between the (numerous) ChPT pa- 
rameters, and the basic QCD coupling and quark 
mass parameters is far from being resolved at 
present. With a more formal (but related) mo- 
tivation, it has been also explored since long ago 
how definite non-perturbative informations may 
be inferred from appropriately modified pertur- 
bation series B, at least in simplified or exactly 



solvable models. In particular, the convergence of 
ordinary perturbation can be systematically im- 
proved by a variational-like procedure, in which 
the separation of the action into "free" and "in- 
teraction" parts is made to depend on a set of 
auxiliary parameters, to be fixed by some opti- 
mization procedure []. 

As a partial attempt to merge some of these 
ideas, we have re-examined |^|| with a new ap- 
proach the above mentioned old problem of gener- 
ating from the basic QCD Lagrangian non-trivial 
values for the quark condensate, pion decay con- 
stant, or dynamical quark masses. The basic 
point is to transform the ordinary perturbative 
expansion, in a s , into an expansion in an arbi- 
trary mass parameter, around a non-trivial (fixed- 
point) solution of the renormalization group evo- 
lution, proportional to the basic scale A. In some 
sense it may be viewed as a systematic, order 
by order, improvement of the original NJL con- 
struction, but with a consistent treatment of the 
renormalization (and directly applied to the QCD 
quark-gluon interactions) . 

2. A crude dynamical mass ansatz 

As a crude first illustration of the mechanism, 
consider the renormalization group (RG) evolu- 



In D = 1 field theories, this optimized perturbation the- 
ory ("delta-expansion") gives a rigorously convergent 
series of approximations, even in strong coupling cases. 



tion of the running mass, 

m(ji ) = m(fj,) exp{~ [ dg 1 ™^ } , (1) 

where (3(g), r ) m (g) drive the running of the cou- 
pling g(fi) and mass m(/i), respectively. Solv- 
ing (|lj) for the "fixed point" boundary condition: 
M = m(M), gives (to first RG order) 



Mi = 



m(/i) 



[l + 26 .g 2 (/i)ln(^)]^ 



(2) 



where bo, 70 are the one- loop RG-coefficients 
(normalization is such that (3(g) — — b g 3 — big 5 — 

■ ■ •> i m (g) = log 2 + 71.9 4 h — )■ 

Although expression (Q) is initially related to 
the "current" mass w(/i) via ([!]), it has the trade- 
marks of a pole mass, thanks to the boundary 
condition M\ = m(M\) ^. This coincidence be- 
tween the pole mass M and the current mass 
m(fi = M), is, however, only an artifact of our 
crude approximation, neglecting at the moment 
the non-logarithmic perturbative corrections Bj. 
Now, the most important property of expres- 
sion (|^) is that it is non-zero in the chiral limit, 
m(n) — > 0. Indeed, (0) identically reads 



Afi(ln(Mi/A))W = rh 



(3) 



where for convenience wc introduced the RG 

_ _ 1 

invariant scale A = p,e 2b os' 2 (at first RG 
order), and the scale-invariant mass rh = 

m(p,)(2b g 2 (fi))~^° ■ ((H) may then be seen as 
a function m(Mi), and requiring its inverse, 
Mi (to), to be defined on the whole physical do- 
main < to < 00, and to match the ordinary 
perturbative asymptotic behavior for to — > 00, 
implies f\ Mi(m — * 0) — > A. It is of course desir- 
able to go beyond the one-loop RG approxima- 
tion, and to take into account as well the non- 
logarithmic corrections, necessary to make con- 
tact with the usual perturbative pole mass 

2 In particular |(3|, ^) is scale invariant (in contrast with 
m(ii)) and gauge-invariant, as the pole mass should be. 

3 Another, a priori possible solution of (H), Mi — * for 
rh — » 0, is rejected because it is only defined for < 
H < (7o/2f>o) 70/ ' 2bo e" 7o/2!,0 A < A, and is therefore not 
compatible with the asymptotic perturbative behavior of 
(§) for m(/i) > A 



Our aim is to obtain a variational "mass gap" 
where the non-trivial chiral limit property of (2|) 
is preserved, while at the same time providing us 
with a systematically (order by order) improv- 
able ansatz, thanks to a particular reorganization 
of the basic perturbative expansion, as will be ex- 
plained in the next section. 

3. Resumming the delta-expansion 

In the present context, a simplest form of 
the so-called delta-expansion Q may be defined 
by formally substituting everywhere in the bare 
QCD Lagrangian: 



too -> m (1 - x); 



9o 



90 (x) 1/2 



(4) 



The parameter x in (Q) just interpolates between 
the free Lagrangian, for x — 0, and the interact- 
ing but massless Lagrangian, for x = 1. In the 
simplest field-theoretical applications, one would 
then use (^) to expand any perturbative expres- 
sion of (mo, <7q) to a given order x q , and try 
to apply some optimization prescription with re- 
spect to the (arbitrary) mass, to . Accordingly, 
the somewhat empirical but most often success- 
ful idea Q] is that the least sensitive region with 
respect to to (entering at any fixed order q) 
should give the best approximation to the exact 
result, which is independent of too. But, in many 
non-trivial field theories, and in particular in the 
present QCD case, before anything the whole pro- 
cedure should be made consistent with renormal- 
ization. As it turns out, the only way to get a 
finite and non-zero result (e.g, M(to — > 0) ^ 0) is 
to resum the a;-series, using an appropriately con- 
structed contour integral transform 0. At first 
RG order, this essentially gives a mass as an inte- 
gral over expression (|^) (with substitution (Q) un- 
derstood). Beyond the one- loop approximation, 
our final mass ansatz reads ^: 



M[(m ) 
A 



2- c m" 
2in 



dv 



F A (v)[C + F(v)] B 
Mi M 2 



F(v) F 2 (v) 



(5) 



4 v in (g) is related to the original expansion parameter x 
as x = 1 — v/q, q being the order of the x-expansion. 



where the contour is around the ] — oo, 0] axis; 

F(v) = \n[m"v] - A InF- (B-C) ln[C + F] , (6) 

withA = 7 1 /(26i),B = 7b/(26o)-7i/(26i),C , = 
&i/(26q); A is the (RG- invariant) scale at two-loop 
order; and finally 



„ _ ( m{n) 
A 



2 C [2b f 



h 2 



1 B 



(7) 



is the scale-invariant, arbitrary (dimensionless) 
"mass" parameter. By construction, F(l) in the 
integrand of (^) resums the leading and next-to- 
leading logarithmic dependence in m (p.) to all or- 
ders |6j . The non-logarithmic perturbative coeffi- 
cients, Mi = (2/3) (7 /26 ) and M 2 , connect [§ 
the pole mass with the running mass m(M). 

Note that it is implicitly always possible to 
choose a renormalization scheme (RS) such that 
bi — ji — for i > 2, since bi, % are then 
RS-dependent. In that sense, eq. (||) resums 
the full RG dependence in ln(m"u). In con- 
trast, the purely perturbative (non-logarithmic) 
information, contained in Mi, Mi, is limited 
by present knowledge to two-loop order. This 
is where the variational principle and optimiza- 
tion play their role, whereby we hope to obtain 
a sensible approximation to the true dynamical 
mass. Observe first that, were we in a simplified 
theory where M\ = M2 = • • • = 0, (0) would 
have a very simple behaviour near its optimum 
(located at to" — > 0), giving a simple pole with 
residue M 2 = (2C)~ C A. Now, in the more re- 
alistic cases, Mi, M2,— cannot be neglected, 
but we can obtain a series of approximants to the 
dynamical mass, by expanding (^) in successive 
powers of m"v, using the standard relation 



— / dv e v v a = -J-^ 
2in J r[-a] 



(8) 



and then looking for optima M^(m" pt ), m" t ^ 0. 
The previous construction is quite general and 
therefore directly applicable to any (asymptoti- 
cally free) model, taking obviously its appropri- 
ate values of the RG coefficients, bi, 7,. The 
ansatz (^) was confronted [^) to the exactly known 
mass gap (TJ for the O(N) Gross-Neveu (GN) 



model, for arbitrary N . The results of different 
optimization prescriptions gave estimates with er- 
rors of 0(5%) or less, depending on N values ||. 
It is important to note that expression ([5]), for 
arbitrary N in the GN model, uses exactly the 
same amount of (perturbative plus RG) infor- 
mation than the one at disposal at present for 
a QCD quark mass: namely, the exact two-loop 
RG-resummed plus perturbative Mi, M2 depen- 
dence. Since our construction essentially relies on 
RG-properties (and analytic continuation), going 
from 2 to 4 dimensions is not expected to cause 
major changes, at least naively. 

4. Hidden singularities of the mass ansatz 

One complication, actually, does occur: as a 
more careful examination of relation (^|) indi- 
cates, there are branch cuts in the v plane, with 
Re[f C Mt] > for the relevant case of rif — 2 or 
3 in QCD. These make the expansion undefined 
when approaching the origin, v = 0, and simply 
indicate the non single-valuedness of (||) below 
those branch points. The origin of those singu- 
larities has some similarity with the renormalon 
ones |Tl"| , as they also appear when extrapolating 
a RG-resummed expression down to an infrared 
scale to" ~ 0. However, a main difference with 
renormalons is that in our construction it is possi- 
ble to move those extra cuts to a safe location, 
Re[i> CMt ] < 0, observing that the actual position of 
those cuts depends, at second order, on the RS, 
via 71 . Performing thus a second-order pertur- 
bative RS change in m(/i), <?(/x), which changes 
7i(M5') to a (singularity-safe) r y 1 , it is then sen- 
sible, in the present context, to invoke a variant of 
the "principle of minimal sensitivity" (PMS) jjj, 
requiring a flat optimum (plateau) of (|5) with re- 
spect to the remnant RS arbitrariness [ 3[ . 
One may perhaps legitimately wonder why the 
ordinary renormalon singularities of the pole 
mass [jfjl do not seem to appear in our construc- 
tion. In fact, the usual renormalon singularities 
always appear as a result of crossing the Landau 
pole 0, which simply reflects an ambiguity from 



Actually, this is an oversimplified picture, valid at one- 
loop RG level only |p^|. However, higher order properties 
of renormalons do not affect, qualitatively, our argument. 



perturbation theory, calling for non-perturbative 
corrections which are typically in the form of 
power corrections In contrast, (||) is such 

that the Landau pole (corresponding to F = in 
our language) is not crossed, but only smoothly 
reached from above, ReF > 0. (Moreover, due 
to the recurrent dependence in F, (|^), implying 
that F(v) ~ to v for m v — > 0, the poles of (||) at 
F = (v = 0) entirely come from the purely per- 
turbative part, i.e. due to M.\,M.2 ^ 0). Note 
that, on more phenomenological grounds, there 
is no strong contradiction with the usual con- 
sequences of the presence of renormalons: while 
the latter indicate, in the pole mass case, an am- 
biguity of 0(A) p^], our construction necessar- 
ily exhibits an arbitrary renormalization scheme 
(RS) dependence, via the above mentioned 71 co- 
efficient, calling for optimization. Practically we 
have obtained: 



in a way similar to the mass case, gives e.g for 



0)~2.97A(2) (9) 
for nf — 2, and a similar result for nf — 3. 



5. Order parameters: F n and (qq) 

The previous dynamical quark mass, although 
it has some meaning as regards DSB in QCD, 
hardly has a direct physical interpretation, e.g. 
as a pole of the S-matrix, due to the confinement. 
In other words, it is not a properly defined or- 
der parameter. It is however possible to apply 
the same construction as the one leading to (|), 
to obtain a determination of the ratios F n / A and 
(qq)(n)/A 3 . The latter gauge-invariant quantities 
are unambiguous order parameters, i.e. F n 
or (qq) ^ imply DSB. The appropriate general- 
ization of (|5|) for Fj. is M 



Fl 2^(m") 2 

P = {2b0) 2z7T 



1 



F 2A-1[ C + F } 



2B 



6„ 1 



dv 
v 

F 



v 2 e' 



F 2 



(10) 



in terms of the same F(v) defined in eq. (g) 
(therefore leading to the same extra cut locations 
as in the mass case), and where S n , and /Sn- 
are fixed by matching the perturbative MS ex- 
pansion, known to 3-loop order 14|. A numerical 
optimization with respect to the RS-dependence, 



F-rr,opt ( m opt 



0) -0.55 A(2) 



(11) 



Concerning (qq), an ansatz similar to ( JlCj ) can 
be derived (with coefficients 5, a, (3 specific to 
(qq) and appropriate changes in the to", F and 
v powers), but for the RG- invariant combina- 
tion m(qq), due to the fact that our construc- 
tion only apply to RG-invariant quantities. To 
extract an estimate of the (scale-dependent) con- 
densate (qq) (/i) is only possible by introducing an 
explicit symmetry-breaking quark mass m exp (i.e. 
mexp 7^ fn), and expanding the m(qq) ansatz to 
first order in m exp . This gives for n/ = 2 |J: 



(qq) l/3 (p = 1 GeV) ~ 0.52 A(2) . 



(12) 



Confronting (|9|) , ( |11| ) and (|12j) gives a fairly small 
value of the quark condensate ^ (and a fairly high 
value of the dynamical mass), as compared to 
other non-perturbative determinations [^5[ . Al- 
though small values of ( qq) are not experimen- 
tally excluded at present m& , it is also clear that 
our relatively crude approximation deserves more 
refinements for more realistic QCD predictions. 

6. Conclusion and discussion 

The variationally improved expansion in arbi- 
trary to", first developed in the GN model ||, 
has been formally extended to the QCD case. It 
gives non-trivial relationships between A and the 
dynamical masses and order parameters, F v and 
(QQ). 

To make progress, what is certainly restrictive 
is the relatively poor knowledge of the purely per- 
turbative part of the expansion (only known to 
two- loop order in most realistic field theories). 
Accordingly, our final numerical results crucially 
depend on the optimization Q. Apart from a few 
models where the series is known to large orders 
(as in the anharmonic oscillator ||,[l7|], or in the 



6 The smallness of (qq) is however essentially correlated 
with the smallness of the F n /A ratio estimate in our frame- 
work, eq. (0). — 

7 For instance, results for (E h, (J1 0[) are substantially dif- 
ferent |p| in the unoptimizeaMa scheme. 



GN model for ./V — > oc), we can hardly com- 
pare successive orders of this expansion to es- 
timate, even qualitatively, the intrinsic error of 
such a method. Invoking the PMS principle jjj, 
although physically motivated, may artificially 
force the series to converge, with no guarantees 
that it is toward the right result. 
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Questions 

E. de Rafael: Can one simply see, just at first 
RG order approximation already, what type of 
graphs are resummed by this mass ansatz? 
J.L.K: It contains the "bubble" chain (the one- 
loop insertions in the gluon line) but, in addition, 
there is an iteration of those dressed gluon lines 
(the so-called Ladder graphs). 
S. Narison: Your calculation is essentially per- 
turbative. Can you include in it the truly non per- 
turbative contributions, like condensates as they 
appear in the operator product expansion typi- 
cally? 

J.L.K: the aim here is to try to estimate these 
NP quantities from the basic QCD interactions 
only, with of course this peculiar resummation. 
To include OPE-like condensates from the start 
would be a kind of double-counting. 



